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Abstract 

The scattering and bound states of the many-body systems, related to the 
short-range Dyson model, are studied. First, we show that the scattering 
states can be realized as coherent states and the scattering Hamiltonian can 
be connected to a free system. Unlike the closely related Calogero-Moser 
model, only a part of the partial waves acquire energy independent phase 
shifts, after scattering. The cause of the same is traced to the reduction in 
the degeneracies. The bound state Hamiltonian for the full-line problem is also 
studied and the relationship of its Hilbert space with that of the decoupled 
oscillators elucidated. Finally, we analyze the related models on circle and 
construct a part of the excitation spectrum through symmetry arguments. 
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I. INTRODUCTION 



Exactly solvable and quantum integrable many-body systems, with long-range inter- 
actions, are one of the most active fields of current research. The Calogero- Sutherland 
model (CSM)liHl, the Sutherland model (SM)§ and their variants are the most prominent 
examples among such systemJl. These models have found application in various branches of 
physicsil"0, ranging from quantum Hall effectinFH, gauge theorieJli'li, chao Jffl, fractional 
statisticsi'&B etc. The CSM and SM type long-range interactions have also manifested in 
models dealing with pairing interactions^ and phase transitionsil. It is known that, the 
CSM is related to random matrix theoryi'i0iHHi and enables one to capture the universal 
aspects of various physical phenomenal. The Brownian motion model of Dyson connects 
the random matrix theory (RMT) with exactly solvable modelsi^. The role of RMT in the 
description of the level statistics of chaotic systems is well-knownili^. 

Some time back, a short-range Dyson model was introduced to understand the spectral 
statistics of systems, which are "non-universal with a universal trend"ill. It is known that, 
there are dynamical systems which are neither chaotic nor integrable, the so called pseudo- 
integrable systems which exhibit the above mentioned level statisticJ^. Aharanov-Bohm 
billiardsii, three dimensional Anderson model at the metal-insulator transition pointS and 
some polygonal billiardsil fit into the above description. Quite recently, a new class of 
one dimensional, exactly solvable many-body quantum mechanical models on the line, with 
nearest and next-to-nearest neighbour interactions, have been introducedlHiil, which are 
related to this short-range Dyson model. Further, using the symmetrized version of this 
model, it has been shown that there exists an off-diagonal long-range order in the system 
which indicates the presence of different quantum phasesil. Apart from possessing a good 
thermodynamic limit i.e., ImiN^ao^ is finite, these models are more physical in the sense 
that, unlike the CSM and SM, where all particles experience pairwise interaction of identical 
strength, irrespective of their distances, here the interactions are only nearest neighbour and 
next-to-nearest neighbour. These models are exactly solvable, but not integrable. Hence, it 
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is interesting to enquire as to how many features of the integrable CSM type systems are 
retained in the present case. For example, the scattering phenomena is quite interesting 
in the CSM case, since the outgoing waves can be shown to be of the incoming type, with 
momenta k'^ = kN+i-i{i = 1, 2, ■ ■ ■ , A^), where fcj's are the incoming momenta. Remarkably, 
the phase shifts are energy independent, a result ascribable to the scale invariance of the 
inverse square interaction. Since, in the present case also, scale invariance holds and the 
interaction goes to zero as the particle separation increases, it is of deep interest to study the 
scattering phenomena. Similarly, for the bound state problem, CSM can be exactly made 
equivalent to a set decoupled oscillators, via a similarity transformationil. Hence, it is also 
of interest to check the same here to understand the precise differences in the Hilbert space 
structure between integrable and non-integrable Hamiltonians, possessing identical spectra. 
It should be mentioned that, in the present case, the degeneracy is less, since the lack of 
quantum integrability i.e., a desired set of mutually commuting operators having common 
eigenfunctions with the Hamiltonian, reduces the degeneracy. Also, these models, being of 
recent origin, need to be analyzed thoroughly, in order to unravel their properties, as has 
been done for the CSM, SM and their generalizations. 

The present paper is devoted to an investigation of these models, and deals with, both 
the scattering and the bound state problems. It is organized as follows: In Sec. II, we 
study the scattering problem and show that the scattering state is a coherent state. The 
connection between the scattering Hamiltonian and that of the free particles in this model 
is then demonstrated. We then study the scattering phase shift and point out its similarities 
and differences with the Calogero case. In Sec. HI, we analyze the relationship of the bound 
state problem with the decoupled oscillators and find some wavefunctions explicitly in the 
Cartesian basis. This analysis reveals explicitly that the degeneracy of this model is less 
as compared to the CSM. Finally, in Sec. IV, the nearest neighbour and next-to- nearest 
neighbour A^-i and BCn models on the circle are studied; a part of their excitation spectra 
is obtained through symmetry arguments. 
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II. THE SCATTERING PROBLEM 



A. Realization of the scattering states as coherent states 

The scattering Hamiltonian, with nearest and next-to-nearest neighbour, inverse square 
interactions, in the units h = m = 1, is given by, 

1 ^ Q2 ^1 ^ 1 

2 ^ dxf ^ ^ ^ {xi — a;j+i)^ ^ (a^j-i — Xi){xi — Xi+i) 

here XN+i = Xi. The corresponding bound state Hamiltonian, contains an additional oscil- 
lator potential: 

It is interesting to note that, the scattering eigenstates can be constructed as coherent states 
of the bound state eigenf unctions like that of the Calogero caseS. To be precise, we show 
that the polynomial part of the bound-state wavefunctions enter into the construction of 
the scattering states. For this purpose, we first identify an SU{1, 1) algebra containing the 
bound and scattering Hamiltonians as its elements, after suitable similarity transformations: 

^ j=l (y^i i=l y^i ~ ^i+l) 



S-\-H/2)S = -i + Eo^ = To . 

Here, Z = Hill I ~ ^i+i 1^ and S = exp{ — ^ Y^i^f} Z exp{ — |T+} . 
Defining, 
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(3) 



2E^NT_ , (4) 



one can easily check that T± and Tq satisfy the usual SU{1, 1) algebra: 

[T+,T_] = -2To , [To,T±] = ±T± . 
The quadratic Casimir for the above algebra is given by. 



C' = T_T+-To(To + l)=T+T_-To(To-l) . (5) 

By finding a canonical conjugate of T+ElhEl: 

[T+,T_] = 1 , (6) 

one can construct tlie colierent state < x \ m, k >, the eigenstate of T+@, wliicli are nothing 
but tlie scattering states, 

<x\m,k>=U-'Prn{x) = e-'^'''^'-PUx) , (7) 

with 



I N Q2 N ^ 

9 H TfT + II 7 — z si^i - 

i=l ^-^i 1=1 ■^i+l) 



Pm{x)=Q. (8) 



It is known that this equation admits homogeneous solutionJ^l z.e, T^Pm^x) = —[{m + 
Eq) /2]Pm{x). Here m refers to the degree of homogeneity of Pm{x). It is easy to see that, 
U^^Pm{,x) is the eigenstate of T+. Starting from Tj^Pm{x) = 0, one gets, 

U-^T+UU-'Pmix) = , (9) 



T+U-'Pm{x) = -^eU~'PUx) . (10) 

The scattering state is given by, ipsca = ZU^^Pm{x) , since Hsca = —ZT^Z^^, therefore, 
Hsca'4'sca = ^^sca- To find < X \ ui, k > explicitly, we have to determine T_. By choosing 
T_ = T_F(To), Eq. (H) becomes 

[T+, T_F(To)] = F(To)T+T_ - F(To + 1)T_T+ = 1 , 
F{To){C + ToiTo - 1)} - F{To + 1){C + T^{T, + 1)} = 1 , (11) 



yielding. 



F{To) = . + " . (12) 
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Here, a is a parameter to be fixed along with the value of the quadratic Casimir C*, by 
demanding that the above commutator is valid in the eigenspace of Tq. Eq. @ when used 
on Pm{x), yields, a = 1 — {Eq + m)/2. Similarly, 



CPUx) = {T_T+ - To(To + l))P„(x) = CPUx) 



where, C = |(m + -E'o)(l — {m + Eq)/2). One then finds. 



(13) 



To - (m + Eo)/2 
Explicitly, we have. 



^ Pm{x) . (14) 



<x\m,k> = e a'^^^-e ^+P™(a;) 



= e-'='/^(fc/2)-(^«-i+^'")(r)-(^«-i+'") JE„_i+„(A;r)P„(x) , (15) 

where, = X^i^;^ and JEo-i+m{kr) is the Bessel function. Note that, there is an additional 
factor of e~"^+ in the above equation, this has been introduced for calculational convenience 
and does not alter our results, since e~^+Pm = Pm- In order to arrive at the above result, 
we have made use of the following: 



T+ (r''^Prn{x)) = 2n{Eo - 1 + m + n)r2("-i)p„(x) 
and also the identityHl, 

oo n 

^oin + a + iy. 

Note that, the above wavefunction can also be obtained by solving the scattering Hamilto- 
nian explicitly. However, we have chosen this algebraic method, since it will be of subsequent 
use. 



B. Connection to Free Particles 



The fact that as in the Calogero case, the spectrum of the scattering Hamiltonian matches 
with that of the free particles and that, the phase shift, as will be shown later, are energy 
independent, suggests a possible connection of this system with free particle^. We now 
show the same by making use of the algebraic structures already introduced. 

The following generators, 

K+ = H,,a, = H,U/3 = 0) , 

1 ^ 

and 

K^ = —Y^{2xA^^) = K^. (16) 

^ i=\ 

satisfy [i^o, K^{K^)\ = ±iK±{K±), [K^K^), K+{K+)] = 2iKo{Ko). The normalization 
of the generators have been chosen differently for convenience. It can be verified that, the 
following operator, 

maps the Hsca, with interactions, to a interaction-free system, i.e., 

uh,,,u^ = hup = o) = -It.^ ■ (17) 

This result motivates one to analyze, explicitly, the precise correspondence of the respective 
Hilbert spaces of the interacting and non-interacting systems. Care has to be taken to ensure 
that, the mapped states are members of the Hilbert space. This analysis, for the present 
case, as well as for the Calogero model has not been carried out so far. We hope to address 
this problem in the near future. 
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C. Analysis of the Scattering Phase Shift 



As is well known, in order to obtain the scattering phase shifts, one has to analyze the 
asymptotic behaviour of the eigenfunctions of the Hamiltonian, when all particles are far 
apart from each other. In the present case, the most general stationary eigenfunction can 
be written as a superposition of all the m dependent states as, 

oo 

^ = ^ E ECm,, (r/2)-(^°-i+'")js,_i+„(fcr)P^,,(a;) , (18) 

m=0 q 

where, Cm,g's are some coefficients and the index q refers to the number of independent Pm's 
at the level m. 

The asymptotic limit of Eq. (|18|) can be obtained from the asymptotic behaviour of 



the Bessel function JEo-i+m{kr) and from the fact that r "^Pm^g{x) does not change in this 
hmit, since Pm,q^ are homogeneous functions of degree m. One obtains, 

^ ~ ^/'in + Vout , (19) 

where, 

\ m=0 q ) 

and 

= e-^(^°-i)-/^(27rfcr)-^ [zr-^^''-^) £ Cm,q (r/2)-™e-™-/2p^,,(x)) e^^^ . (21) 

\ m=0 q ) 

The crucial differences between the Calogero and the present model arise from this point. 
Denoting the equivalent of Pm,,g(a;)'s in the Calogero case, as Pm,g(a^)'s, we first describe the 
Calogero scattering problem, and then compare the results of this model with the same. 
Though the Pm,g(a^)'s, which are the solutions of the generalized Laplace equationi, have 
not been found explicitly so far, for arbitrary m, Calogero had shown the existence and 
completeness of these homogeneous symmetric polynomials. Hence, by choosing suitable 
values for the coefficients C^.^'s, one could characterize the incoming wave, for the Calogero 
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with ki < = 1, 2, ■ ■ ■ , — 1, as the stationary eigenfunction in the center-of-mass 

frame, i.e., J2i^i = 0. Further, since, the Pm,q^ are symmetric under cychc permutations 
and are homogeneous, one obtains, Pm,q{—Tx) = e~^"^'^ Pm,q{x) . where, T denotes a cychc 
permutation of the particle coordinates. Hence, ^/'out for the Calogero model can be written 

as, 

V m=0 q J 

where, k = —k, and —Txi = —XN+i-i- The action of (— T) takes a given particle ordering 
Xi > Xi+i to —X]\r^i_i > —XM-i, and hence preserves the order i.e., the sector of the 
configuration space. This is the reason, invariance of Pm,g(a^)'s under cyclic permutation is 
enough to compute the phase shifts. Comparison with ■?/'in, yields, 

V'out = ce-*(^»-i)- exp |-2 ^ KTx)^ 

-*(^»-i)-exp|z^A;^+i_,a;,| , (23) 



ce 



where, we have used ki = —ki and the cyclic permutation has been carried on the momentum 
variables. From above, it is clear that, the initial scattering situation characterized by 
the initial momenta, ki{i = 1, 2, ■ ■ ■ , A^), goes over to the final configuration, characterized 
by the final momenta. A;- = k^+i-i and the phase shifts are energy independent. In the 
present model, Pm,q{x)^s are the solutions of Eq. (§), which is symmetric only under cyclic 
permutations. Hence, one needs to check if the same steps as narrated above for Calogero's 
case also applies here. First of all, we should find out if the number of Pm,q^ are same 
here. For the sake of clarity, we first consider the four particle case and concentrate on the 
homogeneous solutions of degree four. Since the monomial symmetric functions, provide a 
linearly independent basis set, we can expand P/^{ 

4 4 4 4 4 

P4 = a'^xf + b'^x^Xj + c'^x'fx'^ + d ^ x'^xjXi + e ^ XiXjXiXp . (24) 

i=l ij^j i<j i¥=j¥=^ i<j<l<p 
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The operation of T+ on P4 gives the following three sets of conditions, 

4 

[2(3 + 4:p)a - 2pb + (3 + 4/3)c - 2pd] ^ = , 

i=l 

4 

[Apa + 2(3 + 4/5)6 - 2/5c + 2(1 - p)d - pe] ^ XiXi+i = , 

j=i 

4 

[6(1 + 2/?)6 + 2(1 -2/5)rf] 5^x^x^+2 = , (25) 

where a, b, c, d and e are unknown constants to be determined from the above equations. 
For the purpose of comparison, in the Calogero case. 



2 -^j^ 9xj (xj — Xj) 



Pi{x) = , (26) 



T+P4(x) 
yields, 

4 

[6(1 + 2/?)a - 3pb + 3(1 + 2/?)c - 3pd] ^ x^ = , 

4 

[4/?a + 2(3 + 2/3)6 - 2pc + 2(1 - p)d - pe] ^ x^x^ = . (27) 

i<j 

Hence, the number of solutions are less in the present model as compared to the Calogero 
case. This reduction in the number of solutions takes place because, the interaction term 
contained in the T_|_ operator, leads to the split of the monomial symmetric functions of 
degree (m — 2), giving additional conditions unlike the Calogero case. Explicit calculations 
for a number of few body examples yields similar resultsil. Since the number of Pm,g's also 
represent degeneracy, the same is less here. Hence, it is clear that, Pm.g's do not form a 
complete set. All the solutions obtained so far are symmetric and belongs to a subset of the 
Calogero case. The reason of the symmetric nature of the polynomials lies in the interaction 
term in T+. In order that the action of the interaction term, PYliLi {di — Sj+i), on 

the polynomial Pm{x) results in a polynomial of degree m — 2, the denominator needs to be 
cancelled. This results in the symmetrization of the polynomial, since the nearest neighbour 
couplings arising from ^ connects each particle with every other member of the set. 
The reduction in the number of symmetric polynomials can also be understood from the 
point of view of integrability. Since Calogero model is fully integrable, there are desired 
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number of operators, commuting with the Hamiltonian, which can be used for connecting 
the members of a given set of Pm,(?'s akin to the angular momentum raising and lowering 
operators in the central force problem. The fact that, the number of Pm,g's are less here 
indicates that, the corresponding commuting constants of motion are less here. This point 
will be further elaborated in the next section. 

As has been mentioned earlier, the completeness of the solutions of the generalized 
Laplace equation in the Calogero model enables one to relate all the partial waves of the 
incoming state with those of the outgoing state, with constant energy independent phase 
shifts. As is clear from Eqs. (|20|) , (^1]) and (p2D , in the present case, the outgoing wavefront 
can be made to look like the incoming wavefront. Hovever, only those partial waves, which 
are solutions of Eq. (P) will acquire energy independent phase shifts e"**^-^""^)'^, the rest 
will be unaffected by the interaction. This difference between the Calogero and the present 
model arises, because of the reduction of the number of Pm.g's here. 



III. THE BOUND STATE PROBLEM 
A. Mapping of the model to decoupled oscillators 

The bound state Hamiltonian is given by, 

1 N T N N 1 N 1 

^ ^ - 2 S ^- " i S " ^'^ - ^' g (^^7^ - S - 

where, di = -£r and x^+i = Xi. It is worth pointing out that, for three particles, this model 
is equivalent to the CSM, since the three body term vanishes in this case. The ground-state 
wavefunction and the energy of this system0'0 are respectively given by, iJjq = G Z and 
Eq = {N/2 + NP), where G and Z have been defined earlier. The first term of Eo is the 
ground state energy of oscillators and the second one comes from the interaction. 

In the following, we make use of a method developed in Ref . ( |39|) to show the equivalence 
of this model to a set of decoupled oscillators. For that purpose, we perform a similarity 
transformation on the Hamiltonian, by its ground-state wavefunction, to yield 
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H' = i^o'HiJo = E ^^9^ + Eo-T+ , (29) 

i 

where, T+ = | Y^iLi ^ + PYliLi (x--x-+i) ~ ^^t+i)- Here, we confine ourselves to a sector 
of tlie configuration space given by Xi > X2 > ■ ■ ■ > xjy^i > xjy. Using the identity, 

i 

it is easy to see that 

j 

^From the above diagonahzed form, it is evident that, the spectrum of H is hke that of 
uncoupled oscillators and is linear in the coupling parameter (3. Explicitly, H can be made 
equivalent to the decoupled oscillators: 

1 TV 1 N 

G e-^+('3=o)/2 H e^+('3=°)/^ Q-' = --Y^d^ + -J2 A + ^0 - iV/2 . (31) 

Making inverse similarity transformations, one can write down the raising and lowering 
operators for iY, akin to the CSM. However, the eigenfunctions of H can be constructed 
straightforwardly by making use of Eq. (pUD; since the eigenfunctions of J2iXidi are homo- 
geneous polynomails of degree n in the particle coordinates, n being any integer. Although, 
the similarity transformation formally maps the Hilbert space of the interacting problem to 
that of the free oscillators, it needs a careful study. The singular terms present in T+ may 
yield states, which are not members of the Hilbert space. Below, we clarify this point. 

B. Eigenfunctions in the Cartesian Basis 

In the following, we present some eigenfunctions for the A^-particle case, computed us- 
ing the power-sum basis, Pi{x) = Y.i^=ix[, i.e., ipi = ipoSi; here. Si = e~^+/^P;, and the 
corresponding energy eigenvalue is Ei = {I + Eq), I being an integer. 

The wavefunction (unnormalized) corresponding to the center-of-mass degree of freedom, 
R = jf Y.f=i Xi, is found to be (we use the notation ipni,n2,---,nN = V'o exp{ — iT+} -P"')' 
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1 1 ^ 

V^n.oA- = V^o exp{--T+}i?"^ = exp{--^92} R^^ . (32) 



This can be cast in the form0, 



^ni,o,o,.. = ciJo n I ' ' (33) 

Y^JV i=l '"'I- 

where, Hmi{xi)^s are the Hermite polynomials, and c is a constant. Similarly, the eigenfunc- 
tion for the radial degree of freedom, r'^ = J2 be obtained from, 

^o,n„o,... = ^0 exp{-^T+} (r^)"^ = ^oe-^^+P^ ■ (34) 

For the sake of clarity, we give below the explicit derivation for V'o,i,o,o for the four particle 
case and then generalize the result for arbitrary number of particles and levels. For four 
particles, 

^0,1,0,0 = V^o e-5^+ Pa = V^o e-5^+(x^ + x^ + x^ + x^) . (35) 

We note that, T+Pg = 4(2/3 + 1) , T^Pa = 0, and hence, V^o,i,o,o = - 2(2/? + 1)). For 

particles, it can be verified that, T+r'^"' = 2n(Eo — 1 + n)r^^"~^-' and this gives 4'o,n2fi,- 
aS, 



n2 



m=0 
L 
-'n2 



m!(n2 — m)! (Pq — 1 + ^)! 



c^oL^r'ir') , (36) 



where, L^^'^lr"^) is the Lagurre polynomial. 

Further, for four particles, S3 = P3 — |(2/5 + l)Pi. However, it can be checked that, 
5*4 = e~^+/^P4, does not terminate as a polynomial and results in a function with negative 
powers of the particle coordinates, which is not normalizable with respect to the ground- 
state wavefunction as a measure. This indicates that, there is less degeneracy in the present 
model, as compared to the symmetrized states of the decoupled oscillators. Finding the other 
wavefunctions explicitly, for an arbitrary number of particles, and also the exact degeneracy 
structure of this model remains an open problem. In the above, we have concentrated in 
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finding the wavef unctions in the Cartesian basis; the interested readers are refererred to 
Ref . ( ^6|) for some wavefunctions in the angular basis. 

Below, we list a few eigenfunctions constructed by using the elementary symmetric func- 
tions (we follow the notations of Ref.( ^) for the symmetric polynomials). 

N 

ei = ^ Xi , 62 = ^i^j 5 63 = ^ XiXjXk, ■ ■ ■ ,eN = Y[xi . (37) 

l<i<N i<i<j<N l<i<j<k<N i=l 

In this case, ^{mO = V'o^imj, = e"^+/^ ni(ei)'"% and the corresponding eigen- 

values are, E^^-y = J^i'i^i + Eq. Some of the B^^.ys for the four particle case are listed 
below: 

^2,0,0,0 = ei - 2, 1,0,0 = 6162 - ^(3 - 4/5)ei, ^3,0,0,0 = ^1 - 6ei, 

^2,1,0,0 = ele2 - (3 - 2/5)e? - 2e2 + (3 - 4/3), ^4,0,0,0 = - Uej + 12, 

^3,1,0,0 = e?e2 - ^(9 - 4/5)e? - 66162 + 6(3 - 2/5)ei, ^5,0,0,0 = ef - 20e? + 60ei, 

-64,1,0,0 = e^es - 2(3 - /3)e^ - 126^62 + 12e2 + 6(9 - 4(3)el - 12(3 - 2/3) . (38) 

At this point, it is worth recollecting the Stanley-Macdonald conjecture^, which states 
that, the coefficients of the interaction parameter (3 are positive integers, when the Jack 
polynomialsii are expressed in terms of the monomial symmetric functions with a suitable 
normalization. This conjecture was later proved by Sahiii. Similar feature appears in the 
case of the Hi- Jack polynomialJ^, which are the polynomial part of the wavefunctions of 
the CSM, but with an exception that the coefficient (3 can also be negative. Remarkably, 
from the above explicit computations of the polynomials, we also find that the coefficients 
of the interacting parameter f3 are integers (both positive and negative), though we have 
used elementary symmetric functions. It will be interesting to check whether the modified 
Stanley-Macdonald conjecture also holds in the present case for particles. 



IV. MODEL ON A CIRCLE 



14 



A. An-i Model 

Recently, Jain et al.@i3 have studied a model with nearest and next to nearest neighbour 
interactions and with periodic boundary conditions as given by 



H 



L2 Y sin llixj-Xj+i)] 



f^'^j2 cot[^(xj_i - Xj)] cot[^(a;j - xj+i)] 

3 



(39) 



with Xj+7v = Xi- They have shown that the ground state energy eigenvalue and the eigen- 
function for this model are given by 



N ^ ^2 

xPo = n[sin -(x, - x,+,)f , Eo = Np^ 

3 ^ 



L2 



(40) 



The purpose of this section is to obtain a part of the excitation spectrum of this model. To 
that end, we substitute 



(41) 



in the eigenvalue equation for the Hamiltonian. It is then easily shown that satisfies the 
equation 



1 



N 



N r 



IT 



cot —la;,- — Xj+i) — cot —{Xj^i — Xj 
L L 



dj + Eo-E\(P = 0. (42) 



Introducing, Zj = exp{2iTTXj/L), Eq. reduces to 



HkP = (e - eo) 



(43) 



where 



TV TV 



^3 + 

^3 ~ %+l 



(44) 



with Dj = Zj-^ and e — = (E — Eq)^. It is worth pointing out that the Eqs. (|43| ) and 
(0) are structurally similar to those in the SM except Zk is replaced by Zj+i in our case. 

It may be noted that Hi commutes with the momentum operator P = '^Yl!i=\^i^.- 
Hence is also an eigenstate of the momentum operator, i.e.. 
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P(t) = K(t). (45) 

Further, if is an eigenstate of Hi and P then 

0' = ^^, G = l\^^iz,, (46) 

is also an eigenstate of H' and P with eigenvalues e — eo + Nq^ + 2qK, and k + A^g respectively. 
Here q is any integer (both positive and negative). Note that the multiplication by G 
implements Galilei boost. 

It may be noted that the Hamiltonian and hence the equation is invariant under 
Zj — *• . Since, Zj = e^*'^^'^/^, hence = e~2*'^^j/^ thereby indicating the presence of left 
and right moving modes with momentum k and -n. Hence it follows that, if one obtains a 
solution with momentum k, then by changing zj zj^ ^ one can get another solution with 
the same energy but with the opposite momentum (-/t). Thus all the excited states with 
nonzero momentum are (at least) doubly degenerate. 

Finally, let us discuss the solutions to the equation. So far we have been able to obtain 
the following four solutions. 

(i) = ei , e - eo = 1 + 2/3, 

{ii) = cn-i , t-ei = N -1 + 213, 

{ill) = eiCAT-i - ^ ^ 2if ^' e - eo = + 2 + 4/3, 

iiv) = ctv , e - eo = . (47) 



Here (j=l,2,...,N) denotes the elementary symmetric functions as defined by Eq. 
(defined in terms of Zj). For example, 62 = 2:1^2 + ••• + ^Af-i^Af and it has N( N - l)/2 number 
of terms. 

As mentioned above, each of these solution is doubly degenerate. For example solutions 
ei and e^-i/eN are degenerate. By taking the linear combinations of these two complex 
solutions, it is easily seen that the two degenerate real solutions are 



N N 



^ COS , = XI ' (48) 



i=l i=l 
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where Ui = 27iXi/L. Similarly, all other doubly excited state solutions can be rewritten as 
two independent real solutions. It would appear from this discussion that all the excited 
states are doubly degenerate. However, this is not so. In particular, consider 

= ^1^-^. (49) 

^ 1 + 2/? ^ ^ 

It is easily shown that it is an exact solution to Eq. (^31) with e — eo = 2 + 4/3 but with 
momentum eigenvalue k = 0. This is a nondegenerate solution as it remains invariant under 
Zi In terms of the trignometric functions it can be rewritten as 

(j) = }2 cos{ui - Uj) + . (50) 

i<j n- 

At first sight it appears somewhat surprising that whereas in the Sutherland model, there 
are so many excited state solutions, in our case one is able to obtain so few solutions. In this 
context it may be noted that whereas the Hamiltonian in the Sutherland case is invariant 
under the full permutation group Sn, in our case for > 3 the Hamiltonian Hi as given 
by Eq. has only cyclic symmetry. If one looks at the solutions to Eq. (0), then one 
finds that the condition of no pole in the /3-dependent term almost forces to 'be invariant 
under Sn- Now out of the various Cj (i=l,2,...,N), the only ones in which the demand of 
cyclic invariance necessarily ensures invariance under full permutation group are precisely 
ei, cn-i, ctv, in terms of which we have obtained the four solutions. 



B. BCn Model 



Recently Auberson et al.l23 have studied a BCn model with nearest and next-to- nearest 
neighbour interactions and with periodic boundary conditions as given by 



N 



^ = -2E^ + /5(/5-i)^E 
^ i=i ^-^i ^ i=\ 



2 N 



1 vr 

— TT'- 
sm L 



i+l. 



+ 



1 vr 
sin^ L 



[Xi + Xj+i) 



1=1 



2 N 



TT TT ' 

cot —[Xi^x - Xi) - cot + ^i) 



cot —{Xi - Xi+i) + cot —{Xi + Xi+i, 



TT- ^ 1 TT TT^ ^ 1 27r 



■L2 Vsin' L 
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Following them, we restrict the coordinates Xi to the sector L > xi > X2 > •••• > > 0. 
As shown by Auberson et al.ll, the ground state eigenfunction is given by 



N 



N 



N 



N 



(3/2 



i=l 



i=l 



i=l 



i=l 



where gi, g2 are related to 7, 71 by 



7 

9i = 2 [7 + 271 - 1] , 92 = 271(71 - 1) . 



The corresponding ground state energy turns out to be 



(52) 



(53) 



2L2 



(7 + 271 + 2/3)^ 



(54) 



By setting one or both of the coupling constants 7, 71 to zero we get the other root systems 



I.e. 



5jv : 7i = , Cn -1 = 0, Dn ■■ 1 = li = 



(55) 



The purpose of this subsection is to obtain a part of the excitation spectrum of the 



BCn, Bn,Cn, Dj^ models. To that end, we substitute 



(56) 



in the eigenvalue equation for the above Hamiltonian where i/jq is as given by Eq. 
is easy to show that in that case satisfies the equation 

i:^.2,i:cot,.,-.,«,(55--^) 



21). It 



d 



N Q 

+2/Ji;cot(0,+»,«)(5^+35^ 

+47iEcot20,— + (E-Eo)-^ 
j=i ^^j ^ 



+ 27 ^ cot Oj 



N 

E 

0, 



_d_ 



where 9j = Introducing, Zj = exp{2i9j), Eq. (|F7| ) reduces to 



(57) 



(58) 
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where 

^1 = E + 7 E 7^ A- + 271 E 4^ A- 

j=i j=i i j=i j 

(Dj + Z},+0 . (59) 

Here Dj = zj-^ while e — €q = (E — Eq)^. It is worth pointing out that the above eqation 
is structurally similar to that of the BCj\f of SM, except zj+i is replaced by z^ in our case. 

Note that apart from the cyclic symmetry, the Hamiltonian and hence the equation is 
also invariant under Zj — ^ zj^. As a consequence, as in the BCn Sutherland modelii'0, it 
turns out that even in our case the polynomial eigenf unctions of Hi with BC^ symmetry 
as given by Eq. (|59|) are symmetric polynomials in {zj + j-) i.e. in cos(27rXj/L). 

Finally, let us discuss the solutions to the equation. So far we have been able to obtain 
only one solution in the BCn case but are able to obtain several solutions in the Dn case 
and also a few in the Bn and Cn cases. 



TV 

+/3E 
1=1 



Zj ~\~ Zj 



i+i 

Zj+1_ 



N 

(D,-D,+o + /?E 

7 = 1 



ZjZj+i + 1 



C. Exact Solution for the BCn Model 

It is easily checked that the exact solution is 

(j) = (f)Bc^ = 01 + a , e - eo = 1 + 7 + 271 + 4/3 , (60) 

where 

0, = y(.. + l), a = (61) 

Pi Zj 1+7 + 271 + 4/3 ^ ^ 

We will see that this solution (and in fact other solutions, if any, in the BCn case) will play 
important roles in our construction of solutions for other systems like Bn,Cn and Dn- 



D. Exact Solutions For the Bn Model (71 = 0) 

Apart from the obvious solution (pizj; (3, 7, 71 = 0) as given by Eq. ( |60D it turns out that 
there are other solutions corresponding to the spinorial representation for the Bn model. 
These are 
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nf=i(v^+^), e+-eo = ^[l + 27 + 4/?]. (62) 



In order to obtain the other solution, we start with the ansatz 



= 0+^+ (63) 

where 0"*" is as given by Eq. ( [6^ ) and consider the equation Hicp = (e — eo)0. It is easily 
seen that in that case satisfies 

Ht{z,; P, 7, 71 = 0)^+ = (e - 6+)^^+ , (64) 
where is as given by Eq. (^) while 

^ z-1 

Ht = H.izf, 7, 71 = 0) + -^D, 

j=i + ^ 

= i/i(z,;/?,7-l,7i = l). (65) 

We thus see that ip'^ essentially satisfies the same equation as satisfied by the BC^ Hamil- 
tonian but with the value of 7 and 71 shifted to 7 — 1 and 1 respectively. Thus it follows 
that once we obtain solutions of the BCn problem, all these will give us new solutions of 
the spinorial type for the Bj^i model as given by Eqs. (|6^) to (|65|). Unfortunately, so far we 
have been able to obtain only one solution in the BC^ case as given by Eqs. (^) and (^TJ). 
Using that solution, it then follows that the new spinorial solution for the S^v model is as 
given by Eq. ( |OBD with energy 

N 

e-eo = 2 + 7 + 4/3+-[l + 27 + 4/3], (66) 
while ip^ = (pBcA^j] /5, 7 - 1, 7i = 1) with (pBCi, being given by Eqs. (ID and (|l|). 

E. Exact Solutions For the Model (7 = 71 = 0) 

Apart from the above solutions (|60|) , (|6^) and (|63D (with 7 = 71 = 0), we have found 
several other solutions in the Djsi case. This is related to the fact that unlike Bj^, there are 
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two distinct classes of spinor representations for D^. Besides, there are also some additional 
solutions in this case. It may be noted that in this case the Hamiltonian Hi acting on is 

Hiizf, (3) = Hiizj- /3, 7 = 0, 71 = 0) . (67) 

The first new solution that we have is given by 

nf=i(v^ - ^) , e- - eo = ^[1 + 4/5] . (68) 



Note that the two solutions ( |62D (with 7 = 0) and (|68| ) which correspond to the two different 
spinorial representations are degenerate in energy. 

In order to obtain the other solution, we start with the ansatz 

^ = <t)-il}- (69) 

where is as given by Eq. (|68| ) and consider the equation Hicj) = (e — eo)*/). It is easily 



seen that in that case ip satisfies 

//r(%;/3,7 = 0,71 = 0)V'" = (e - e-)r , (70) 
where e~ is as given by Eq. (^) while 

^ Z- + 1 

Hi = Hiiz,; /3, 7 = 0, 71 = 0) + ^ ^^T^^' 

j=i 

= ^ri(2:,;/5,7 = l,7i = 0). (71) 

We thus see that ip' essentially satisfies the same equation as satisfied by the Bj^i Hamiltonian 
but with the value of 7 being fixed at 1. Using the solution for the case as given by Eqs. 
(|60|) and (|6lD (with 7 = 1, 71 = 0) it then follows that the new spinorial solution for the 
model is as given by Eq. (|69D with energy 

e-eo = 2 + 4/3 + -[l + 4/3], (72) 

while ip~ = (pBCNi^jj /^5 7 = 1)71 = 0) with (PbCn being given by Eqs. (^) and (pl|). Notice 
that this solution is degenerate in energy with the solution (|6^) (with 7 = 0). 
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In addition, we find that the product of the two "spinorial solutions" is also a solution 
of the Djq model, i.e. 

= 0+0- = nf=,(z, - -) , e+- - eo = N[l + 2/5] . (73) 

In order to obtain another solution, as above we start with the ansatz 

= 0+0-^+- (74) 

where 0+, 0~ are as given by Eqs. (|62|) and ( |67|) respectively and consider the equation 
Hi(f) = (e — eo)0. It is easily seen that in that case ip^ satisfies 

Ht-{zy, /?, 7 = 0, 71 = 0)^+- = (e - e+-)^+- , (75) 

where is as given by Eq. ([73|) while 

^ + 1 

= /fi(z,; /5, 7 = 0, 71 = 0) + 2 5] ^^T^^' 

= i/i(^,;/?,7 = 0,71 = 1). (76) 

We thus see that ■0"'"" essentially satisfies the same equation as satisfied by the Cn Hamil- 
tonian but with the value of 71 being fixed at 1. Using the solution for the BCn case as 
given by Eqs. (|60D and (|6lD (with 7 = 0,71 = 1), it then follows that the new solution for 
the Dm model is as given by Eq. ( |74|) with energy 

e - eo = 3 + 4/5 + Ar[l + 2/3] , (77) 
while = (f)BCNi^j'^ 7 = 0, 7i = 1) with 0BCiv being given by eqs. (|0D and (|6l|). 



F. Exact Solution For the Cn=4 Model (7 = 0) 

So far we have discussed all the solutions which are valid for any N. In addition, in the 
special case of = 4, we have been able to obtain a solution in the CAr=4 (and two solutions 
in the 0^=4) case. The solution for the 0^='^ case is given by 
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a31 



A 



(^1 + -)(^2 + -)(^3 + -) + CP] + 501 

Zl Z2 Z3 



where 



e - eo = 3 + 671 + 8/5 , B 



1 + 271 + 2p 

Here, by CP. one means cyclic permutations and 0i is as given by Eq. (|60|). 



(78) 



(79) 



G. Exact Solutions For the Dn=4 Model (7 = 7^ = 0) 



Clearly, one obvious solution in the Dm=a case is obtained from the solution (|78| ) by 
putting 7i = 0. The other solution is obtained by making use of the ansatz as given by Eq. 
(0), i.e. let 



a31. 



(80) 



with ((r"^ being given by Eq. (^). Using Eqs. (|7^ to ( [TGD it then follows that this is a 
solution for the -DAr=4 model with 

8A/3 



e-eo = 13 + 16/3, B 



3 + 2/5 



(81) 



V. CONCLUSIONS 

In conclusion, we have carried out a systematic study of the many-body Hamiltonian, 
related to the short range Dyson model. The scattering state of this model is obtained 
and is shown to be a coherent state. Akin to the CSM, the connection of the scattering 
Hamiltonian to a free system is established. Unlike the Calogero model, analysis of the 
scattering process for the present model reveals that, only a part of the partial waves acquire 
energy independent phase shifts. We then showed the mapping of the bound system to 
decoupled oscillators and by explicitly computing some of the bound-state eigenf unctions, 
we find that, the present model has less degeneracy as compared to the Calogero-Sutherland 

23 



model. Finally, we have studied the ^Iat-i, BCn, Bn, Cn and Dn models on a circle and 
obtained a part of their excitation spectrum. A number of open problems, like finding the 
Pm,q{xys explicitly for the N-body problem, characterizing the degeneracy structure, finding 
the complete eigenspectra and conserved quantities for these type of models still remains to 
be tackled. We hope to come back to some of these issues in future. 

M.E, N.G and P.K.P would like to thank Prof. V. Srinivasan and Prof. S. Chaturvedi 
for useful discussions. 
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